Infrared catastrophe and tunneling into strongly correlated electron systems: 

Perturbative x-ray edge limit 
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The tunneling density of states exhibits anomalies (cusps, algebraic suppressions, and pseudogaps) 
at the Fermi energy in a wide variety of low-dimensional and strongly correlated electron systems. 
We argue that in many cases these spectral anomalies are caused by an infrared catastrophe in the 
screening response to the sudden introduction of a new electron into the system during a tunneling 
event. A nonperturbative functional-integral method is introduced to account for this effect, making 
use of methods developed for the x-ray edge singularity problem. The formalism is applicable to 
lattice or continuum models of any dimensionality, with or without translational invariance. An 
approximate version of the technique is applied to the ID electron gas and the 2D Hall fluid, 
yielding qualitatively correct results. 

PACS numbers: 71.10.Pm, 71.27.+a, 73.43.Jn 



I. TUNNELING AND THE INFRARED 
CATASTROPHE 

The tunneling density of states (DOS) is known to ex- 
hibit spectral anomalies such as cusps, algebraic suppres- 
sions, and pseudogaps at the Fermi energy in a wide va- 
riety of low-dimensional and strongly correlated electron 
systems, including 

(i) all ID metalsj 1 ! 2 ! 3 ! 4 ! 5 

(ii) the 2D diffusive metali 6 i 7 i 8 i 9 i 10 i n 

(iii) the 2D Hall fluidj 12 i 13 i 14 i 15 . 16 i 17 i 18 i 19 i 20 i 21 

(iv) the edge of the confined Hall 

fl mc l 22.23.24.25.26.27.28.29.30.31.32.33.34.35.36.37.38.39.40.41 



In this paper we propose a unified explanation for these 
anomalies, develop a nonperturbative functional-integral 
formalism for calculating the electron propagator that 
captures the essential low-energy physics, and apply a 
simplified approximate version of the method to cases (i) 
and (iii), yielding qualitatively correct results. 

We claim that the physical origin of the DOS anoma- 
lies in the above systems is the infrared catastrophe of 
the host electron gas caused by the sudden perturba- 
tion produced by an electron when added to the sys- 
tem during a tunneling event. This infrared catastro- 
phe, a singular screening response of a conductor to a 
localized potential turned on abruptly in time caused by 
the large number of electron- hole pairs made available by 
the presence of a sharp Fermi surface, is known to be re- 
sponsible for the singular x-ray optical and photoemission 
spectra of metals j 5 ^ 5 ^* 5 ^ Anderson's related orthogonal- 
ity catastrophe; 57 ! 58 and the Kondo effecti 59 ! 60 To un- 
derstand the connection to tunneling, imagine the tun- 
neling electron being replaced by a negatively charged, 
distinguishable particle with mass M. In the M — > oo 
limit, the potential produced by the tunneling particle 
is identical — up to a sign — to the abruptly turned-on 



hole potential of the x-ray edge problem, and an infrared 
catastrophe of the host electrons would be expected. 

Tunneling of a real, finite-mass electron is different be- 
cause it recoils, softening the potential produced. How- 
ever, in the four cases listed above, there is some dy- 
namical effect that suppresses recoil and produces a po- 
tential similar to that of the infinite-mass limit: In case 
(i) the dimensionality of the system makes charge relax- 
ation slow; In case (ii), the disorder suppresses charge 
relaxation; In case (iii), the Lorentz force keeps the in- 
jected charge localized; and case (iv) is a combination 
of cases (i) and (iii). It is therefore reasonable to expect 
remnants of the infinite-mass behavior. The real electron 
is also an indistinguishable fermion, unable to tunnel into 
the occupied states below e-p. 
,42,43.fiM.46?»ai : h/|js 50 ,pj.§£<geds as follows: We use a 
^rmbbard-StratonovTcn^ransformation to obtain an exact 
functional-integral representation for the imaginary-time 
Green's function 

G(rf<7f,riC7i,ro) = -(T^ CTf (r f ,r o )^ ffl (ri,0)) H , (1) 

and identify a "dangerous" scalar field configuration 
<Axr(r, t) that causes an infrared catastrophe. Assum- 
ing To > and an electron-electron interaction potential 
U(r), " 



xr (r, t)=U[t- R(t)] e(r)O(T - r) 



(2) 



where R(r) = rj + ( rt rc| ri )r is the straight-line trajectory 
connecting ri to r t with velocity (rt —t{)/tq. For the case 
of the tunneling DOS at point r , we have r; = r f = r 
and 



<Mr,r) = tf(r - r )e(T)e(T - r), 



(3) 



which is the potential that would be produced by the 
added particle in Q if it had an infinite mass. Fluctua- 
tions about </> xr account for the recoil of the finite-mass 
tunneling electron. 

In this paper we introduce an approximation where 
the fluctuations about </> X r are entirely neglected, which 
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we shall refer to as the x-ray edge limit. We expect the 
x-ray edge limit to give qualitatively correct results in the 
systems (i) through (iv) listed above; the effect of fluctu- 
ations will be addressed in future work. The x-ray edge 
limit of our formalism can be implemented exactly, with- 
out further approximation, only for a few specific models. 
Here we will instead carry out an approximate but more 
generally applicable analysis of the x-ray edge limit by 
resumming a divergent perturbation series caused by the 
infrared catastrophe. We apply this "perturbative" x-ray 
edge method to the ID electron gas and 2D spin-polarized 
Hall fluid, both having a short-range interaction of the 
form 



The Euclidean propagator can be written in the 
interaction representation with respect to Hq as 



(TV^ t .(r f ,r )y; (Ti (r i ,0) e -/o i ^^)) 

(Te-IodrV(r) )o 

AM Dfj,[4>] g(rfer f ,ri(7i,To|0), (7) 



U{r) = XS(r). 



(4) 



The exact implementation of the x-ray edge limit for 
these two cases will be presented elsewhere. 



II. INTERACTING PROPAGATOR 

We consider a D-dimensional interacting electron sys- 
tem, possibly in an external magnetic field. The grand- 
canonical Hamiltonian is 



n_ 

2m 



v (r) - ^0 



G(rf£rf,ri£ri,T " 



where 

Here J\f= (Texp(— dr V))q 1 is a To-independent con- 
stant, and 

ff(rf<7j,ri<Ti,To|0) 
= - <TV^ f (r f ,To)^ i (r i ,0)e^o' 3rf ^ rfD '-^ fa ( r ^) (9) 

is a correlation function describing noninteracting elec- 
trons in the presence of a purely imaginary scalar poten- 
tial — i(p(r, t). In the following we always assume that 
to > 0. 

Next we deform the contour of the functional integral 
by making the substitution <f> — > i(f> XI + 4>, leading to£i 



Mr) 



G(r f CTf , r iCTi , T ) = A^.f <^-V xr 
+ i^y'd £5 rd D r>t( r )^t /(r / )t/(r _ r / ) ^ /(r ' ) ^ (r)j xy > D M [0] e -^^ _1 ^(?(r f( 7 f ,r i( T i ,To|# xr + 0). (10) 



where II = p+ - A, and where vo(r) is any single-particle 
potential energy, which may include a periodic lattice 
potential or disorder or both. Apart from an additive 
constant we can write H as Hq + V , where 



H = j2j dD ^i( r ) 



2m 



-I- v(t) - fl 



Mr) (5) 



and 



V=\ d D rd D r' 8n(r)U(r-r')Sn(r'). 



(6) 



Hq is the Hamiltonian in the Hartree approximation. The 
single-particle potential v(r) includes the Hartree inter- 
action with the self-consistent density no(r), 

v(r) = v (r) + I d D r' U{r-r') n (r'), 



and the chemical potential has been shifted by —U(0)/2. 
In a translationally invariant system the equilibrium den- 
sity is unaffected by interactions, but in a disordered or 
inhomogeneous system it will be necessary to distinguish 
between the approximate Hartree and the exact density 
distributions. V is written in terms of the density fluc- 
tuation Sn(r) = V4( r )Vv(r) - n (r). 



This is an exact representation for the interacting Green's 
function, where <f> now describes the fluctuations of the 
Hubbard-Stratonovich field about <A xr . 



III. DOS IN THE X-RAY EDGE LIMIT 

As explained above, our method involves identifying 
a certain field configuration tfi xr that would be the po- 
tential produced by the tunneling particle if it had an 
infinite mass. Apart from a sign change, this potential is 
the same as that caused by a localized hole in the valence 
band of an optically excited metal. As is well known from 
work on the x-ray edge singularities, such fields cause an 
infrared catastrophe in the screening response. Fluctua- 
tions about 4> XT account for the recoil of a real, finite-mass 
tunneling electron. 

In this section we introduce an approximation where 
these fluctuations are entirely neglected, which we shall 
refer to as the extreme x-ray edge limit. The approxi- 
mation can itself be implemented in two different ways, 
perturbatively in the sense of Mahan£4 and exactly in 
the sense of Nozieres and De Dominicis^ In Sec. IIVI 
we apply the perturbative x-ray edge method to the ID 
electron gas, including spin, and to the 2D spin-polarized 
Hall fluid. The x-ray edge method will be implemented 
exactly for these and other models in future publications. 
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In the x-ray edge limit, we ignore fluctuations about Eqs. (|llfl and (|12fl define the interacting propagator in 
in ^(rftTf , rjCi, T O |i0 xr + (ft), approximating it by the x-ray edge limit. The local tunneling DOS at position 



g(rfCTf, i-jCT;, r \i(ft xr )r& Then we obtain, from Eq. (fTTTfl . 

G(r f CT f , i-iCTi, t ) = NgirfOi , riCi, T \i<ft XI ), (11) 
where, according to (J5J, 



r is obtained by setting r; = rf = r and ct; = <jf = ctq, 
and summing over ctq. 



^(rflTf^ilTijTol^xr) 

-(T^ t (r f ,r )^ i (r i ,0)e^^ 5 ") . 



(12) 



J 



A. Infrared catastrophe and divergence of perturbation theory 

To establish the connection between Eqs. (|llfl and l|12l) and the x-ray edge problem, we first calculate the tunneling 
DOS at ro by evaluating (|12f> perturbatively in <ft XI , 



5(roCTo,r CTo,To|# X r) = G (r a-o, r o°o, to) + J d n dn </>xr(ri,n) G (r cro, rier , t - n) G (ria , r <To, n) 

Go(r cro,r 1 cr ,To - n) GoC'lOo, ^oo, n - r 2 ) G (r 2 cr , r cr , t 2 ) 
0(C), (13) 



+ / d D rid D r2dr 1 dT2<^a.(r 1 ,ri)<^ a .(r 2 ,T2) 
- i G (rocro,rocro,To)no(ri,r 2 ,Ti - t 2 



Here Go(rcr, rV, t) = — (T^ CT (r, 7")^v(r',0)}o is the mean-field Green's function associated with i?0: an d 

n (r,r',r) = -(TSn(v, r) Sn(r', 0)) 

is the density-density correlation function, given by Co(rcr, r'cr, r) Go(r'cr, rcr, — r). Now we use (0, and assume 
the short-range interaction Q. Then Eq. (|13[) reduces to 

g(r ao, "* (To, To|#xr) = G (r cr, r cr, r ) + A / dn G (r cr , r fT , to - n) G (r cro, r oo, n) 

Jo 

G (rocro,r cro,To - n) G (r <7o, roao, n - r 2 ) G (r a , r a , r 2 ) - § Go(r o- , r ao, r ) n (r , r , n - r 2 

(14) 



A / dr% dr2 
Jo 

0(A 3 )- 



We evaluate (1141) for the ID electron gas and the 2D Hall fluid in the large to£f limit using the low-energy propagators 
of Appendix^] e F is the Fermi energy. In the ID electron gas case this leads to 



g(r ao,r ao,T Q \i(ft XI ) = G (r a , r a , t )\ 1 - 2AAT ln(r e F ) + A 2 7V ( 2 



2 In (r e F ) - 2 ln(r e F ) + — tq6f 



. (15) 



where we have used the asymptotic results given in Appendix [B] and have kept only the corrections through order A 2 
that diverge in the large ro limit. Here No is the noninteracting DOS per spin component at e F S These divergences 
are caused by the infrared catastrophe and the associated breakdown of perturbation theory at low energies. Similarly, 
for the Hall fluid we find 



g(r o o- o ,r o cro,To|i0xr) = Go(r cr , r f7o, t )< 1 - A 



1 — v\ 1 , o / 1 — V 
TO + T.^ 



2ir£ 2 



2tt£ 2 



(16) 



where v is the filling factor and t is the magnetic length. The divergence in this case is stronger because of the infinite 
compressibility of the fractionally filled Landau level at mean field level. 



B. Perturbation series resummation 

A logarithmically divergent perturbation series simi- 
lar to (|15f) occurs in the x-ray edge problem, where it 



is known that qualitatively correct results are obtained 
by reorganizing the series into a second-order cumulant 
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expansion. Here we will carry out such a resummation 
for both and JTSJ. 

In the electron gas case this leads to 



g(r ao,r a ,T Q \i(j> XY ) w G (r Q a Q , r Q a , r ) 



/ 1 



where 



\ T oeF 



a = 2XN + 2A 2 7V 2 . 



,f A 2 Af 2 r e F 



(17) 



(18) 



The exponential factor in i|17|) . after analytic continua- 
tion, produces a negative energy shift. However, this shift 
depends sensitively on the short-time regularization and 
is not reliably calculated with our method, a is equal to 
the x-ray absorption/emisssion edge exponent 

2(<p/n) + 2^/tt) 2 

of Nozieres and De Dominicis^ (including spin) for a 
repulsive potential, with ip = arctan(7rAiVo) the phase 
shift at cf, when expanded to order A 2 . 
In the Hall case with < v < 1, 



where 



g(r a ,r a ,T \i^ x 



7 



v — 1 
2ttP 



,7(i/-1)t 



2^ 2 



(19) 



(20) 



is an interaction strength with dimensions of energy. This 
time the infrared catastrophe causes a positive energy 
shift and no transient relaxation. The energy shift in this 
case is predominantly determined by long-time dynamics 
and is physically meaningful. 



generic algebraic DOS of the Tomonaga-Luttinger liquid 
phase is at least consistent with our assertion that the 
algebraic DOS in ID metals is caused by the infrared 
catastrophe. 

We note that the x-ray edge approximation actually 
predicts a power-law DOS in clean 2D and 3D electron 
systems in zero field as well. However, in those cases 
there is no reason to expect the x-ray edge limit to be 
relevant: In those cases the electron recoil and hence 
fluctuations about </> xr are large. 



B. 2D Hall fluid 

For the 2D Hall fluid with < v < 1 we obtain 

N(e) = const xS(e-[l- v]-y) . (22) 
The lowest Landau is moved up in energy by an amount 



A 



2ir£ 2 ' 



(23) 



The DOS (|22|) is also qualitatively correct, in the fol- 
lowing sense: The actual DOS in this system is thought 
to be a broadened peak at an energy of about e 2 / k£, with 
k the dielectric constant of the host semiconductor, pro- 
ducing a pseudogap at ep. Of course the magnetic field 
dependence of the peak positions are different, but the 
physical system has a screened Coulomb interaction, and 
and here we obtain a hard "gap" of size (|23|l at ep. We 
speculate that the cumulant expansion gives the exact 
x-ray edge result for this model, but that fluctuations 
about <fr XI will broaden the peak in l|22|) . 



DISCUSSION 



IV. APPLICATIONS OF THE PERTURB ATIVE 
X-RAY EDGE LIMIT 

The tunneling DOS is obtained by analytic continua- 
tion. We take the zero of energy to be the mean-field 
Fermi energy. 

A. ID electron gas 

For the ID electron gas case we find a power-law 

N(e) = const xe Q , (21) 

where the exponent a is given in (|18|) . In l|21|l we have 
neglected the energy shift appearing in (|17|) . 

The power-law DOS we obtain is qualitatively correct, 
although the value of the exponent certainly is not. The 
exponent will be modified by carrying out the x-ray edge 
limit exactly, and also presumably by including fluctua- 
tions about xr . However, the fact that we recover the 



The results and 1)21(1 and H22J1 are consistent with our 
claim that the DOS anomalies in the ID electron gas and 
the 2D spin-polarized Hall fluid are caused by an infrared 
catastrophe, similar to that responsible for the singular 
x-ray spectra of metals, the orthogonality catastrophe, 
and the Kondo effect. (The models we considered are 
of course tractable by existing specialized methods.) In 
addition to providing a common explanation for a vari- 
ety of tunneling anomalies, our method may provide a 
means of calculating the DOS in other strongly corre- 
lated and low-dimensional systems, such as at the edge 
of the sharply confined Hall fluid investigated experimen- 
tally by Grayson et al.jQ Chang et al.J* and by Hilke et 
ai.pa* where existing theoretical methods fail. 

Anderson has gone even further and proposed that this 
same infrared catastrophe causes a complete breakdown 
of Fermi liquid theory in 2D electron systems, such as the 
Hubbard model, in zero fieldi£i££*£& We are at present 
unable to address this question with the methods de- 
scribed here, which assume that the effects of fluctuations 
about xr are small. 
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Our method can be applied, in principle, to realistic 
tunneling geometries, with or without disorder, and with 
Coulomb interaction. 



APPENDIX A: LOW-ENERGY 
NONINTERACTING PROPAGATOR 

For a noninteracting electron system with single- 
particle states </>a(r) and spectrum e Q , the imaginary- 
time Green's function defined in Eq. Q is given by 

G (ra, rV, r) = <W ]T ^(r)0;(r')e-( e °-^ r 

a 

x (jn F (e a - M) - 1] + « F (e Q - m) 6(-t)) , (Al) 

where np(e) = (e' 3 ' 5 + l)" 1 is the Fermi distribution func- 
tion. In this appendix we shall evaluate the diagonal 
components Go (rcr, rcr, r) for two models in the large |t|, 
asymptotic limit. 



ture with filling factor v. In the gauge A = Bxe y , 

Mr) = c nk e~ ik \ e-Wlt-U? Hn ( x/e H), (A5) 
where c n k = (2 n ri}.TT^£L) 2. Here t = y/hc/eB is the 
magnetic length and L is the system size in the y di- 
rection. The spectrum is e n = huj c (n +5), with lu c = 
eB/m*c the cyclotron frequency (to* is the band mass). 
In the lowest Landau level < v < 1, 

4> k = (tt^L)^ e- lky e -WI^)\ ( A6 ) 

At long times r ^ uj~ x we obtain 

v-Q{t) 



Go(r, r, r) 



2irP 



(A7) 



APPENDIX B: ASYMPTOTIC EVALUATION OF 
TIME INTEGRALS 



1. ID electron gas 

The first model is a translationally invariant electron 
gas at zero temperature in ID (the derivation we give is 
actually valid for any dimension D). In the limit of large 
cpM we obtain 



G (r<7, re, r) 



No 

5 

r 



(A2) 



where iVo is the noninteracting DOS per spin component 
at the Fermi energy It will be necessary to regu- 

larize the unphysical short-time behavior in Eq. i|A2|l . 
The precise method of regularization will not affect our 
final results of interest, such as exponents, which are de- 
termined by the long-time behavior. We will take the 
regularized asymptotic propagator to be 



G (rcr, rcr,r) 



Re 



iVo 



t + i/e F 



(A3) 



When possible, we will let the short-time cutoff in 
approach zero, in which case (|A3|) simplifies to 



G (rcr, r<7,r) 



p iVo 

T 



(A4) 



where P denotes the principal part. The results (|A3|I and 
(|A4|) are valid for any spatial dimension D; the only D 
dependence appears in the value of Nq for a given ep. 



2. 2D Hall fluid 

The second model we consider is a 2D spin-polarized 
electron gas in the lowest Landau level at zero tempera- 



Here we note the asymptotics used to obtain (I15fl 
1 \ „ / 1 \ 2 



dr Re — Re . 

u V T 0- r + V e F/ V r + V e F 



1 



Re 



x Re 



drdr' Re. , ... , 

\t - T + i/e F J \T-T' + i/e F 

1 



T 



1 



ln(r e F ) 



t' + i/e F 



— ln 2 (T e F ), 

TO 



and 



drdr' 



Re 



1 



t — t' + i/eF 



- T e F - 21n(r e F ). 



In these expressions we have retained all terms, including 
subdominant contributions, that diverge in the to — > 00 
limit. 
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